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Abstract 

In a thermal field theory, the cumulants of the momentum distribution can be extracted 
from the dependence of the Euclidean path integral on a shift in the fields built into the 
temporal boundary condition. When combined with the Ward identities associated with 
the invariance of the theory under the Poincare group, thermodynamic potentials such 
as the entropy or the pressure can be directly inferred from the response of the system 
to the shift. Crucially the argument holds, up to harmless finite-size and discretization 
effects, even if translational and rotational invariance are broken to a discrete subgroup 
of finite shifts and rotations such as in a lattice box. The formulas are thus applicable at 
finite lattice spacing and volume provided the derivatives are replaced by their discrete 
counterpart, and no additive or multiplicative ultraviolet-divergent renormalizations are 
needed to take the continuum limit. In this paper we present a complete derivation of 
the relevant formulas in the scalar field theory, where several technical complications are 
avoided with respect to gauge theories. As a by-product we obtain a recursion relation 
among the cumulants of the momentum distribution, and formulae for finite- volume 
corrections to several well-known thermodynamic identities. 



1 Introduction 



Thermal field theory is the theoretical tool for computing properties of matter at high 
temperatures and densities from first principles [HE]- It allows one, for instance, to 
determine the equation of state of Quantum Chromodynamics, which in turn is an es- 
sential ingredient to understand the properties of matter created in heavy ion collisions, 
and to model the behavior of hot matter in the early universe (for recent reviews see 
Ref. [llll]). Obtaining first-principles predictions from a thermal field theory is often 
challenging since it describes an infinite number of degrees of freedom subject to both 
quantum and thermal fluctuations. Even though there are several established methods 
to compute the thermal properties of field theories [SHS] , new theoretical concepts and 
more efficient computational techniques are still needed in many contexts particularly 
when weak-coupling methods are inapplicable. 

In a recent Letter we related the generating function of the cumulants of the to- 
tal momentum distribution to a path integral with properly chosen shifted boundary 
conditions in the compact direction normalized to the ordinary thermal one [9]. By 
exploiting the Ward Identities (Wis) associated with the space-time invariances of the 
continuum theory, the cumulants can be related in a simple manner to thermodynamic 
potentials. In a relativistic theory at zero chemical potential, for instance, the variance 
of the momentum probability distribution measures the entropy of the system. Thanks 
to a recursion relation among the cumulants, the kurtosis is related to the specific heat. 
These results suggested a new way to determine the equation of state of a thermal field 
theory [9]. 

The aim of this paper is to present a complete and self-contained derivation of the 
formulas introduced in Ref. [9| for the scalar field theory, where several technical com- 
plications in their derivation are avoided with respect to gauge theories jlQ]. A crucial 
role is played by the symmetry-constrained path integrals |llU12j and by the continuum 
Wis associated with the relativistic invariance of the theory. The latter allow us to 
derive the recursive relation among the cumulants, to generalize well-known thermo- 
dynamic relations to finite- volume systems, and to exclude additive or multiplicative 
ultraviolet-divergent renormalizations of the cumulants. The formulas are applicable, 
up to harmless finite-size and discretization effects, at finite volume and lattice spacing, 
where ratios of path integrals can be determined by ab initio Monte Carlo computa- 
tions [9] . As a result the entropy density, the pressure and the specific heat of a thermal 
field theory can be obtained by studying the response of the system to the bare shift 
parameter. 

After a section on the basic properties of the scalar theory, the symmetry-constrained 
path integral and its relation with the cumulant generator is introduced in section [3l 
The relevant Wis are derived in sections H] to [H and the main results and conclusions 
of the paper are reported in section [71 Several technical details are given in the four 
appendices. 
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2 Preliminaries and basic notation 



We are interested in the thermal scalar theory defined via the Euclidean path integral 
formalism in an infinite volume as well as in a finite box of volum^ V = with the 
field (j) satisfying periodic boundary conditions. The time extent is set to Lq = 1/T, 
where T is the temperature of the system. Some basic definitions and Wis, associated 
with the invariance of the continuum theory under the Poincare group, are reviewed in 
this section. Other properties of the theory, not directly relevant for the subject of this 
paper, can be found in various textbooks, see for instance Refs |13yi4j. 
The partition function of the theory is defined as usual as 

Z = j D(j)e-^ , (2.1) 

where the action S = f d^x C is defined by the Lagrangian density 

C = \ {d,m^.<t>) + V{<i>) , y(0) = \ mV' + ^ • (2.2) 
The equations of motion are given by 

||^-n0|(x)Oi...O„,^ =f^(Oi...O^...O„), (2.3) 

where Oj is a generic local field inserted at the point a:*, and 0[ is its variation with 
respect to the fundamental field and its derivative at the coordinate value x. 

2.1 Translational invariance 

The theory is invariant under space-time translations, i.e. under the direct product of 
four continuum Abelian groups of global transformations 

x'^, = x^,-e^, , ^\x'^) = 0(x^) . (2.4) 

The associated Wis can be derived in the usual way by promoting the symmetry to a 
local one, i.e. — )• £^{x), and studying the variation of the functional integral under 
the local transformations 

S(f){x) = £f,{x) df,(j){x) , 5{d,,cl){x)} = d^{e^{x)df,(j){x)} . (2.5) 

By considering different functions Ef^iz) and properly chosen fields Oj, many interesting 
non-trivial relations can be derived. For ei,{z) = e^S^'^^z — x) one obtains 

n 

eu (d^^T^uix) Oi . . . 0„) = - ^ (Oi . . . 5f O, ...On) , (2.6) 

1=1 



^Throughout the paper the linear dimension in the spatial direction k will be indicated by Lk- 
^Sometimes the argument is shown explicitly to clarify the meaning of some formulas. Summation 

over repeated indices is understood unless explicitly specified. No summation over k is understood in 

sections 13] to [S] and everywhere for the field Tkk. 
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where 6^0i is the variation of the field Oj under the transformation (|2.5p . and the field 

V = (9^<A)(5.(A) - (2-7) 

is the energy- momentum tensor of the theory symmetric under the exchange jj, -v^ v. 
The Wis ()2.6p . and therefore the consequences discussed in the following, applies to 
connected correlation functions as well. When all operators Oi are localized far away 
from X, the classical conservation identities 

(a^r^,(x)Oi...o„) = (2.8) 



are recovered. If we integrate Eq. (12. 6p over a bounded region R which contains in 
its interior the points x^, . . . , x™" while the fields Om+i ■ ■ - On are localized outside, the 
Gauss theorem leads to 

„ m 

eu / da^{x) {T^u{x)Oi . . . 0„) = - V (Oi . . . 6,0, ...On) , (2.9) 

where (5eOj is the variation of the field Oi under the global transformation associated 
with e^, and the integration measure da^{x) points along the outward normal to the 
surface dR. In this case, the very same Wis could have been obtained by considering 
directly the limit where ei,{x) goes to a constant in i?, and it is null outside. 

The discussion above applies also to the theory defined in a finite box with periodic 
boundary conditions, provided the (5-function is replaced by its periodic generalization 
(J^P-* . In the following we will use the same notation for both of them since the precise 
meaning will be clear from the context. 

2.2 Time evolution of operators 

From the Wis in Eq. ()2.9p we can derive the time-evolution of a generic field. If we 
choose = SuqEq, and we integrate over a thick time-slice R with the field Oi being 
an arbitrary functional of (j) and d^(j) inserted into it while the operators O2 • • • On are 
localized outside, then 

(^Oi(x) = eo5oOi(x) (2.10) 

and therefore 

do{0i{x^)02...0n) = - [ dao(x)(roo(x)Oi(xi)02...0n ) . (2.11) 

JdR 

This is the Euclidean version of the time-evolution of the generic field Oi(x^). The field 

Tooixo) = j cZ^xToo(x) (2.12) 
is therefore (minus) the Hamiltonian of the system, and we can write 

(Too) = J-lnZ, (2.13) 
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or more generally 

(Too(Lo)0)e = ^(0), (2.14) 

where O is a generic operator which does not depend explicitly on the time coordinate 
and it is at a physical distance from the time-slice Lq. 

2.3 Translational invariance of correlators 

If all fields Oi,...,0 n are inside the integration region R, then the transformation (|2.5p 
implies 

60i{x) = e^d^0i{x) (2.15) 

for arbitrary functionals of (p and d^(f). In the limit where the bounded region R goes to 
infinity or, for instance, in a finite volume with periodic boundary conditions, the l.h.s. 
of Eq. (|2.9p can be neglectecH. Correlation functions enjoy translational invariance, i.e. 
they satisfy 

n 

^ {Oi{x') . . . Oi{x') . . . 0„(x'^)) = . (2.16) 

i=l 

As an example of transformation of a composite field, it is interesting to study the 
variation of the energy-momentum tensor under local transformations. If we rewrite 
the tensor (|2.7p as 

T^u = ^ RfiuaP {da(t>){dl3<t>) " S/iuVic/)) , R^vaP = ^^aKfi + S^/S^ua " S^^Sa/^ , (2.17) 

its variation under a local space-time translation (|2.5p can be written as 

6T^i,{x) = ep{x) dpTpu{x) + R^uap {Tppix) + 5pi3C{x)} da£p{x) , (2.18) 
and Eq. (j2.15p is satisfied if £p{x) is constant. 

2.4 Invariance under 4-dimensional rotations 

In the Euclidean the invariance under the homogeneous Lorentz group is replaced by 
the symmetry under 5*0(4) rotations. An infinitesimal transformation reads 

x'„= Xa+tJapXis (2.19) 

with ojap = —ijJjia and, as in the previous section, we can derive the associated Wis by 
studying the variation of the functional integral under an infinitesimal local transfor- 
mation. The latter can be written as 

x'a = Xa- Eaix) , £aix) = -UJa/six) X/s , (2.20) 

^In Wis associated to space-time translational invariance no boundary terms of this type are expected 
to contribute at infinity in the quantum theory, at variance of what happen in those associated with 
spontaneously broken symmetries. 
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and the fields transform accordingly to Eq. (j2.5p . If we choose Ua/s (x) = Wap [z — x], 
we obtain 

n 

wap {d^K^,^p{x) Oi . . . 0„) = -2 ^ (Oi . . . . . . 0„) , (2.21) 

i=l 

where 

K^-aP = XaT^p - Xi3Tna (2.22) 

is an antisymmetric tensor for a o /3. The Wis (j2.2ip are just combinations of those 
in Eqs. (|2.6p since, thanks to the symmetry of T^^, d^K^.^p is a linear combination 
of the components of d^T^^p with field-independent coefficients . All Wis associated 
with 4-dimensional rotations can thus be derived from those discussed in the previous 
section. It is, however, instructive to consider properly chosen integrated Wis. They 
can be obtained by integrating over a bounded region R which contains in its interior 
the points x^, . . . , x*" while the fields Om+i ■ ■ ■ On are localized outside. This leads to 

„ m 

w^fi / d^x {d^K^,^p{x) Oi...On) = -2j2{Oi--- S^Oi ...On) (2.23) 
•^^ i=i 

where 5«,0j is the variation of Oj under the global transformations associated with Wa^. 
In particular it is interesting to consider an infinitesimal boost, that in the Euclidean 
is obtained by choosing coQj = vj and all other components null, an integration region 
R which is a thick-time slice between the two hyper-planes with xq = i/q i: t, and a 
correlation function with a field Tofc(yo) inserted inside R while the operators Oi . . . On 
are localized outside. The variation of space-time components of the energy-momentum 
tensor under a global boost reads 

STokiy) = Vj^yo djTokiy) - yj 9o?ofe(y) + Tjkiy) - SjkToo{y)^ , (2.24) 
and the integrated Wis can thus be written as 

/ daf,{x) {Kf,.oj{x) Tokiy) Oi...On)= (2.25) 
JdR 

-(|yo9jTofc(y) - yjdoTokiy) + Tjkiy) - SjkToo{y)^ d . . . 0„) . 

By using Eq. (j2.8p . and for correlation functions for which boundary terms can be 
neglected when integrating by parts in d^y, such as connected correlations functions 
with operators Oi localized in space, the Wis read 

[ da^ix) {K^.oj{x)Tok{yo) d • • • On)c = Sjk (Too (2/0) Oi • • • On)c • (2.26) 

JdR 

Being valid for every string of external localized operators, i.e. for a generic state, this 
is the Euclidean version of the commutation relation of the momentum with the charges 
associated with the boosts. Finally it is interesting to notice that the number of fields 
Tfj^i, inserted in the correlation functions entering the two sides of these Wis is different. 
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3 Momentum distribution from shifted boundaries 



In this section we show how the correlation functions of the total momentum can be 
extracted in the Euclidean path integral formalism by generalizing the periodic temporal 
boundary condition to shifted boundary conditions. The relative contribution to the 
partition function of the states with momentum p is {Lq dependence suppressed) 

R{p) ^ Tr{e-^«Hp(p)| 
V Tr{e-^oH} ' 

where the trace is over all the states of the Hilbert space, P^p^ is the projector onto 
those states with total momentum p, and H is the Hamiltonian of the theory. If we 
introduce the partition function 

Z{z) = Tr{e-^o^ e'^'} , (3.2) 

in which states of momentum p are weighted by a phase e*P'^, and we use the standard 
group theory machinery (see Ref. [12] for a detailed discussion on this point) it easy to 
show that 

i?(p) = ^ / d^ze-'^-'Z{z), (3.3) 

where Z = Z(0) is the ordinary thermal partition function. The generating function 
K(z) of the cumulants of the momentum distribution is defined as usual as 

e-''^^^ = ^^e'P-^Rip), (3.4) 
p 

and the cumulants are given by 

t.v ^ ( 1 '|ni+W2+n3+l ^^""^ "^^"^ ^^"^ ^('^) 

{2ni,2n2,2n3} ^ 5zf ^ Sz^"^ ^zf ^ V z=0 ' ^ ' 

where they have been normalized so as to have a finite limit when V —?■ oo. The shifted 
partition function Z{z) can be expressed as a Euclidean path integral with the field 
satisfying the boundary conditions 

(/>(Lo,x) = (/.(0,x + z) (3.6) 



in the compact direction. From Eqs. (|3.3p and (j3.4p . the generating function can thus 
be written as the ratio of partition functions 

e-^(^) = ^ , (3.7) 

i.e. two path integrals with the same action but different boundary conditions, and 
the cumulants can be obtained by deriving it with respect to the shift parameter z an 
appropriate number of times. Since the cumulants are connected correlation functions 
of the total momentum charge of the theory, they are ultraviolet finite as they stand 
(see also section [5] and appendix |A]). Therefore K(z) and the momentum distribution 
R{p) are finite as well, and they do not need any ultra-violet renormalization. 
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3.1 Extension to the lattice 



When defined on a lattice, the theory is invariant under a discrete subgroup of trans- 
lations and rotations only, the momenta are quantized, and the continuum Wis are 
broken by discretization effects. Generic lattice definitions of the energy-momentum 
tensor, as well as the corresponding charges, require ultraviolet renormalization. It is 
still possible, however, to factorize the Hilbert space of the lattice theory in sectors with 
definite conserved total momentum. The formula for the lattice projector P^p) can be 
easily obtained, and its explicit form can be found in Ref. |12j . Since only physical 
states contribute to the symmetry constrained path integrals in Eq. (|3.ip . the lattice 
momentum distribution R{p) is expected to converge to the continuum universal one 
without the need for any ultra-violet renormalization. The definition of the cumulants in 
Eq. (|3.5p is thus applicable at finite lattice spacing, provided the derivatives are replaced 
with their discrete counterpart, and no additive or multiplicative ultraviolet-divergent 
renormalization is needed for taking the continuum limit. 



4 Two-point correlators of T. 



MM 



The two-point correlation functions of the diagonal components of the energy-momentum 
tensor satisfy Wis which are non-trivial and interesting in the thermal theory at finite 
volume. By choosing T^kiu) interpolating operator, = d^QeQ and by using transla- 
tional invariance and parity, the WI ()2.6p can be written as 



&'o{{Tm{xo)Tkk{y))+5{xo-yQ) (Too + /:)} = . (4.1) 

Analogously by taking = 6^k^k, Tqq{z) as interpolating operator, and by using trans- 
lational invariance and "/c-parity" 

d]:{{fkk{wk) Tooiz)) + 5{wk - Zk) {Tkk + /:)}= , (4.2) 



where 

T^u{wk)= I \\\dwp T^y{w) . (4.3) 

By integrating both of them and by taking the difference, we obtain (xq ^ yo ; ^^fc 7^ ^k) 

Lo {Too{xo) Tkk{y)) - Lk {Tkk{wk) Tqq{z)) = (Too) - {Tkk) ■ (4.4) 

If we subtract from both terms on the l.h.s the disconnected piece, we arrive at 

Lo {Toq{xo) Tkk{y))c - Lk {Tkk{wk) Tqq{z))c = (Too) - {Tkk) ■ (4.5) 

Non-trivial properties of Eqs. (j4.4p - (|4.5p are that all operators are localized at non- 
zero physical distance from each other, that the number of fields T^^ inserted in the 
correlation functions entering the two sides of each equation is different, and that the 
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additive renormalization counter-terms proportional to ai and 02 in Eq. (|X3]) do not 
contribute. These equations are thus vahd for the renormahzed energy-momentum ten- 
sor as well, or conversely they imply that the overall renormalization constant satisfies 
Z = 1 (see appendix [A]). In regularizations that break translational invariance, they 
can be imposed to compute non-perturbatively the overall renormalization constant of 
the diagonal component of the field T^i^. 



4.1 Thermodynamic limit and finite-size effects 



In the thermodynamic limit the second term on the l.h.s. of Eq. (j4.5p vanishes, and the 
infinite volume WI reads 



Lo {Too{xo)Tkk{y))c = (Too) - {Tkk) 



(4.6) 



The very same WI can be obtained in a slightly more elegant form. By integrating a 
combination of two of the Wis in Eq. (j2.6p , and by using translational invariance and 
parities we obtain 



(fx (xo d^T^o{x) Tkk{y) - xk d^T^k{x) Tqq{z)) = (Too) - {Tkk) 



(4.7) 



from which it is easy to derive Eq. ()4.6p . The latter has a straightforward thermody- 
namic interpretation. If we remember that in the Euclidean 



(Too) 



{Tkk) = P 



(4.8) 



where e and p are the energy density and the pressure in the thermodynamic limit 
respectively, Eqs. (|2.14p and (|4.6p imply the well-known thermodynamic relation 



^ dp 



By defining the entropy density of the finite-volume system as usual 

Tlnzl , 



VdT 



(4.9) 



(4.10) 



and by using Eq. (|2.14p and the analogous one in the k-direction (see Eq. IB.4P ). the 
extensivity of the free energy and Eq. (14. 9p lead to 



Ts = e + p 



d_ 
dT 



P 



(4.11) 



The leading finite-size corrections in these thermodynamic potentials were calculated in 
Ref. [15j for a generic theory with a mass gap M in the screening spectrum. For the 
sum of the energy and pressure they read 



{Tkk) - {Too) - (e + p) 



uMT 
' 27rL 



M + 3T- 



BM 



dT 



(4.12) 
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where the factor v stands for the multiphcity of the hghtest screening state, and the 
dots stand for terms which vanish with a larger exponent. The WI (|4.5p generahzes 
Eq. (I4.9P to finite-volume, and it shows that the violations in the latter are due to the 
finite-volume dependence of the energy density, i.e. 



(Tfcfc) - (Too) - T 



d{Tkk) 



d{Too) 



dT dLk ' 

By inserting the finite-size corrections computed in Ref. [15] we obtain 



-t — — {J-kk) + {J-oo) - 



dM 



1 + 



1 



dT ^ ' ^ 27r \dT J [' ' ML 

and analogously for the violation of first equation (j4.1ip 



(Tkk) + (Too) 



3 + 3ML + M'^L^ 



(4.13) 



(4.14) 



(4.15) 



Finite-volume effects in the thermodynamic relations (j4.9p and (j4.1ip are exponentially 
small in ML, and the leading finite-size corrections are known functions of the lightest 
screening mass, its temperature derivative and the multiplicity of the associated state. 



5 Two-point correlators of Tq^ 

In a finite box a generic boost transformation rotates periodic fields to non-periodic 
ones. Those transformations are incompatible with the boundary conditions of the 
finite- volume theory. The integrated Wis associated with the 50(4) rotations must 
be modified by finite-size contributions which vanishes in the thermodynamic limit. 
On the other hand the finite-volume theory has an energy-momentum tensor which is 
locally conserved, and the integrated Wis associated with a generic infinitesimal rotation 
UajB can be constructed starting from Eqs. p.6p and by multiplying both sides by the 
proper factors. By choosing = 6,^k^k, Tofc(y) as interpolating operator, and by using 
translational invariance and parity, the WI (j2.6p can be written as 

dE{{Tok(.xo)Tok{y)) - S{xo - yo) (T^fc + /:)} = . (5.1) 

Analogously by taking = 5uoeo, TQk{z) as interpolating operator, and by using trans- 
lational invariance and parity 

d^{{fok{wk)Tok{z)) - 6{wk - zk) (Too + C)] = . (5.2) 
By integrating both of them and by taking the difference we obtair0 {xq ^ yo ^ z^) 
Lq (Tofc(xo) Tok{y)) - Lk{fok{wk) TQk{z)) = (Too) - (Tkk) ■ (5.3) 

^This relation is explicitly verified for the free theory in appendix [Cl 
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This WI can be derived directly from Eq. (j2.25p by paying attention to the fact that 
boundary terms in the integration by parts cannot be neglected. The comments below 
Eq. (I4.6P apply also in this case. In particular, in regularizations that break translational 
invariance, these equations can be imposed to compute non-perturbatively the overall 
renormalization constant of the off-diagonal components of the field T^j^. These Wis are 
reminiscent of those associated with non-singlet chiral symmetry in QCD, which turn 
out to be so relevant to renormalize the axial current when chiral symmetry is explicitly 
broken by the regularization |161ll7j (see also Ref. [18j for a review). From Eqs. (j4.5p 
and ()5.3p we obtain 

{Tokixo)Tok{y)) - {TooixQ)Tkkiy'))c = j^\^{Tok{wk)Tokiz)) - {fkkiw'k)Tooiz'))c^ , 

(5.4) 

where in all correlators the operators must be inserted at a physical distance (xq ^ i/q, 
etc.). It is interesting to notice that this relation can be derived directly from the 
conservation of the energy-momentum tensor in Eq. (12. 8p without explicit reference to 
the transformation properties of the interpolating operators, see Ref. [9j and below. 

5.1 Thermodynamic limit and finite-size effects 

In the infinite-volume limit the second term on the l.h.s. of Eq. (j5.3p vanishes, and the 
relation reads 

Lo (rofe(xo) Tokiy)) = (Too) - {Tkk) ■ (5.5) 

Its thermodynamic interpretation is straightforward. If we remember that in the Eu- 
clidean the momentum operator maps to pk —iTq^, then (xq / yo) 

{To3(.xo)To3{y)) = -^'{0,0,2} ■ (5.6) 

The Eq. (15. 5p can then be written as 

k{o,o,2}=T{e + p) = Th . (5.7) 

The WI (15. 3p generalizes this equation to finite- volume. The second term on the left- 
hand side of (15. Sp vanishes exponentially in the lowest screening level corresponding 
to a state with a non-zero momentum in the time direction. Whenever the lightest 
screening state has vanishing momentum in the time direction, which is expected to 
be the (generic) case, the leading finite-volume corrections to the second cumulant are 
those from the r.h.s. of the Eq. ()5.3p . They are exponentially small in ML, and their 
explicit form is given by the r.h.s. of Eq. (|4.12p multiplied by T. 

6 Recursion relation for 2n-point correlators of Tq^ 

In the thermodynamic limit the Wis (|2.6p imply a recursion relation among correlators of 
Tofe inserted at a physical distance. It can be derived by repeatedly using Eq. ()2.6p with 
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different strings of interpolating operators. This can be concisely shown by introducing 
a new action 

Sj = S- [ d^xJ{xo)Tok{x) (6.1) 



which differs from the standard one by a term which couples an external source J{xq), 
constant over the time-slices, with the momentum field in direction k. The correspond- 
ing path-integral is defined as 

-Sj 



Z[J] = j Dcj) e~^J , (6.2) 

and the expectation value of a generic operator {0)j is defined as usual. It turns out 
to be useful to introduce the operators ^.n] defined as 

where F(J) is a generic functional of the external source and the Xq are all different. 
As usual when applied to F(J) = lnZ[J] it gives the connected correlation functions of 
n momentum operators Tofc(xo) inserted at a physical distance. By choosing = S^k^k 
(no summation over A;), Tqo as interpolating operator, where 

/ Yl dxp\ T^^ix) , (6.4) 

and by using translational invariance in tim^, the WI (|2.6p gives 

dfD[,^^^,2,.^{%o{x')fok{x^))j,c = 9f Z)[,3..,2„](ffcfe(:E2)foo(x^))j,c . (6.5) 

Analogously by taking = 6^oeo, Tqj. as interpolating operator, and thanks to trans- 
lational invariance and the symmetry of Tq^, the WI (|2.6|) leads also to 

dfD^^^^^^2.^{%oix')fok{x^))j,c = 9f Z?[,3..,2„j(Tofc(xi)Tofe(x2))j,c . (6.6) 
By putting together last two Wis we arrive at 

5f D[,3..,2„]{(Tofc(:r2)To,(xi))j,e - {fkk{x^)%o{x'))j,c} = . (6.7) 

Since the argument of the partial derivative is constant in x^, we can integrate by 
keeping all insertions at a physical distance (xg all different) and obtain 

Dixl.x^^^]{{Tok{xl)Tok{xl))j^c - {Tooixl)Tkk{xl))j,c'j = 

Ll i^[^3..^2„]|(rofe(xi)rofc(x2))j,c- {Tooix^)Tkkix'^))j,c} , (6.8) 



5 

is the standard one 



When an operator D^^i ^2^] is applied, the action entering the definition of the correlation functions 
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which generahzes Eq ()5.4p . The Eqs. (jB.SP and (jB.lip give 

'^[xl.xl^']{Lo{Tkk{xl))j - Lk-^lnZ[J] - J dwo J{wo) {Toki'Wo))j ,^ = (6.9) 

where again the Xq are all different, and (a generalization of) Eq. ()2.14p leads to 

— — d — 

D[xl.^2rr]{Too{xl)Tkk{xl)).ic = ^ L>[^3.^2„](rfcfe(x^))j . (6.10) 

Finally, by putting together the last three equations, we can write (all insertions at a 
physical distance from each other) 

{TokixD . ■■ToUxln)c = (2n - 1) ^ {i- {Tok{xl) ■ ■■Tokixln)c} + 
Ll{{fok{x')fok{x^)Tok{xl) . ..Tok{xln)c - {foo{x')fkk{x^)Tok{xl) . . . Tofc(x2")),} 

{Tokixl)...Tok{xln)c} . (6.11) 
6.1 Thermodynamic limit and finite-size effects 

The last three terms in Eq. (j6.1ip are finite-size effects which vanish in the infinite 
volume limit: the first two because the distance — x'^\ can be arbitrarily large, the 
third one due to the expected volume dependence of the correlation function. In the 
thermodynamic limit we thus arrive at the recursive relation 

{TokixD . ..Tokixln)c = (2n - 1) ^ {i- {Tok{xl) . . . ro,(xg"))e} , (6.12) 
with a straightforward thermodynamic interpretation 

^{0,0,2n} = (2"- - l)7'^^{r/i:{o,0,2n-2}} • (6.13) 

Cumulants with non-trivial indices in the other two spatial directions are related to 
those in Eq. (|6.13p by cubic symmetry. For the total momentum, Eq. (|6.12p is the 
analog of the well-known recursive relation among the connected correlation functions 
of the energy, see Eq. (|2.14p with O being a string of Tqc's. The Eq. (j6.13p is checked 
explicitly in appendix [D] for the free theory. It is interesting to notice that it leads to a 
straightforward physical interpretation of the fourth cumulant, i.e. 

k{om = 3r^c, +9T\e + p), (6.14) 

where Cy is the specific heat of the system in the thermodynamic limit. Finite- volume 
corrections to this identity can be computed by starting from the Wis ()B.ip and by 



u 
'dT,, 
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following a procedure analogous to the one which lead to Eq. ()5.3p . For the fourth 
cumulant the relevant WI read^ 

{(Too)-(Tfc,)} (6.15) 

—3 |Lo {TqqTqq)c — Lk {TkkTkk)c^ + 6 {TQkTokToo)c - L\ {T^kTQkTkk) ■ 
To compute the leading finite-size effects, it is convenient to rewrite this equation as 

-^^0 {TokTokTokTok)c = ^{-^0^^ ^^'dLk ~ ^}{^"^°°'* ~ C^kk)^ 

^ ^ r Q d ^ 

{TokTokTokTok)c — 6Lfc|l + ^k-gj^ — ^o-gj-j {TokTok) ■ (6.16) 

The second line of the right-hand side vanishes exponentially in the lowest screening 
level corresponding to a state with non-zero momentum in the time direction. We 
therefore turn our attention to the first line, from which we compute the leading finite- 
size corrections to the fourth cumulant by using again the results in Ref. [15] 

{TokTokTokTok)c - 3Tc, " + = "^{^ ^ (e'^^^MT^ + 



g-ML ^rj. ,^ _ 2^^) /g ^ g^^^ ^ M^L^) \ + ... (6.17) 

oT > 

where the dots stand for terms that vanish exponentially faster. 
7 Main results and conclusions 

By comparing the WI (|5.7|) and the Eq. (j3.5|) with = 1 and n\=n2 = 0, the entropy 
density in the thermodynamic limit can be written as 

Thanks to Eq. (j4.1ip . the pressure can be computed by integrating s in the temperature, 
and the ambiguity left due to the integration constant is consistent with the fact that 
p is defined up to an arbitrary additive renormalization constant. From the relation 
()6.14p . the specific heat of the system can be written as 



-ML 



= lim ^ 

y-5-oo V 



2 1 



+ 7^- 



lnZ({0,0,z}) . (7.2) 

2 = 



The last two equations make clear that the response of the partition function to the shift 
z is governed by basic thermodynamic properties of the system, and that the potentials 



In this and next equation the arguments Xg and x]. are omitted. Operators are inserted always at 
a physical distance. 
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entering the equation of state of the thermal theory can be extracted by rather simple 
formulas. It is worth stressing that Eqs. (j7.ip and (j7.2p are the result of a judicious 
combination of the Wis associated with invariance of the theory under the Poincare 
group only. The relevant Wis relate correlation functions of conserved charges, and this 
is why they are free from ultraviolet subtractions and renormalizations. The formulas 
()7.ip and (j7.2p are indeed valid for a wider class of thermal field theories than the scalar 
one [9j. For gauge theories, however, their derivation involves additional complications 
due to the non-commutation of translations and gauge transformations, and it will be 
presented in a forthcoming paper |10) . 

Crucially the Eqs. (17. ip and (j7.2p remain valid in a lattice box, up to exponentially 
suppressed finite-size effects and harmless discretization errors, provided the derivatives 
are replaced with their discrete counterpart. In a finite volume some of the relevant 
Wis generalize well-known thermodynamic relations, and they are the basic ingredient 
to relate the finite-size effects in the cumulants to those in simple quantities such as 
the energy density and pressure which are known [15]. With respect to the analogous 
observables computed in the standard methods [SHH], the entropy density and the spe- 
cific heat computed from the equations above do require neither a vacuum subtraction 
nor an ultraviolet renormalization constant to be fixed. Moreover an improvement of 
the action automatically leads to a corresponding improvement in the thermodynamic 
quantities. Conversely the shifted boundary conditions adopted in this paper and the 
Wis (14. Sp and ()5.3p . which are valid for the renormalized energy- momentum tensor as 
well, can be combined to desig n a non-perturbative renormalization procedure for 
on the lattice. 

The generating function K(z) may turn out to be an interesting thermodynamic 
quantity in itself. For a scale-invariant theory, for instance, one might have expected 
that K{z) can be an arbitrary function of T\z\. However, the recurrence relation ()6.13p 
and the scale invariance fix its functional form unambiguously to be 



no matter what the coupling of the theory is. This is particularly relevant to those 
strongly coupled theories that can be treated with the AdS/CFT correspondence. We 
note that for z purely imaginary, v = iTz can be interpreted as the macroscopic velocity 
of the thermal system. The double poles in z appearing in Eq. (j7.3p corresponding to 
|v| = 1 therefore do not come as a surprise. 
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A Renormalization pattern of 



The renormalization pattern of the energy-momentum tensor was studied in great detail 
in perturbation theory |19H23j . Here we assume to work in a regularization which 
preserves translational invariance, such as the dimensional one [ 19ti22j . A more general 
analysis, which is needed for the lattice theory can be found in [23]. Since the tp"^ 
theory is not asymptotically free, it is difficult to extend the results reviewed here non- 
perturbatively. For the scope of this appendix, however, the standard power counting 
argument is sufficient. 

The bare energy-momentum tensor is a dimension-4 field, symmetric under the 
exchange ^ ^ u, even under the transformation — > —cj), and which satisfies the 
conservation equation (12. 8p . The renormalized tensor can thus be written as 

%^ = Z {T^^ + ai [5^^,U - ^^^^;\ (P^ + aa • (A.l) 

Since by construction the last two fields on the r.h.s. have zero divergence, T^i, satisfies 
the conservation equations (j2.8p as well. If we choose the Oi to be a combination of 
elementary fields, the integrated Wis (|2.9p can be written as 

„ n 

% \ da^ix) {%,{x)di . . . 0„) = - V (Oi . . . 5 A ...On) , (A.2) 
2 JdR ^ \ / 

where Oi are the corresponding renormalized fields, and the contributions proportional 
to ai and 02 vanish because they are surface integrals of vectors with null divergence. 
Since the correlation functions entering Eq. ()A.2p are finite by construction, then also 
2 must be finite. If we choose Z = 1, the renormalized energy-momentum tensor reads 

f^,y = T^u + ai [SfiuO - d^du] 0^ + a2 S^,y (A. 3) 

and it satisfies Wis of the same form as those of the bare field, i.e. Eqs. (|2.8p and (|2.9p . 
The last two fields on the r.h.s. of Eq. ()A.3P do not contribute to the charges Tofc(xo), 
since the integration by parts of the term proportional to ai gives no boundary contri- 
butions. The physical momentum fields are therefore defined by their bare expressions, 
and no renormalization ambiguity is left. The field proportional to 02 contributes to 
the charge Too(a;o) by an additive constant term 02^^, which is ultraviolet divergent 
and proportional to the volume of the system. It contributes to the vacuum expectation 
value of the energy, but it cancels in connected correlation functions of Tqo with other 
fields. An analogous analysis applies for the fields T}.y{xk) {v 7^ k) and T/^kixk). They 
are the momenta and the energy if direction k is interpreted as the "time" direction. In 
particular the divergent coefficient 02 is the very same for Too(a;o) s^nd T^kixk). 

The second term on the r.h.s. of Eq. (1A.3P (as well as the third one) contributes 
to the trace T^^ of the energy momentum-tensor. In the quantum theory the renormal- 
ization procedure indeed breaks the conformal invariance of the classical theory also in 
the massless limit (trace anomaly). 
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B Ward identities in presence of external sources 



In this appendix we generalize some of the results of section [2] to the theory defined by 
the action Sj given in Eq. (16. ip . The starting point are the the Wis 

{^d^T^uix) + J {xo)dyTQk{x) - do [J{xo)iT^k{x) + 6^k^x))] 

n 

-J{xo)dk[T,o{x) + 6Mx)]] Oi . . . On) J = - ^ (Oi . . . 6^0i 

i=l 

which extend those in Eq. (12. 6p . 

B.l Evolution and translational invariance in A;-direction 

From the previous Wis we can derive the "k-time" evolution of a generic field. If we 
choose e,y = 5^k^k (no summation over k) in Eq. (|B.ip . and we integrate over a thick k- 
slice R with the field Oi being inserted into it while the operators O2 • • ■ On are localized 
outside, we obtain 

dk {Oi{x') O2 . . . On) J = - [ dak{x) {Tkkix) Oi{x') O2 . . . On )j ■ (B.2) 

JdR 

This is the Euclidean version of the fc-time evolution of the generic field Oi{x^), and 
therefore the field T^k is the Hamiltonian of the system associated with the fc-time. If 
R covers the full space, we obtain 

n 

^ df (Oi(xi) . . . 0,(x*) . . . Onixn)j = , (B.3) 

i=l 

i.e. a generic correlation function is translational invariant in the k direction even in 
presence of the external source. This is expected since the source field is constant 
in Xk, and it does not break translational invariance in this direction. By following 
an analogous procedure, it is straightforward to show that translational invariance is 
preserved in the other two spatial directions too. Since the field Tkk is the Hamiltonian 
in direction k, we can write 

{fkk{xk))j = In Z[J], (B.4) 
dLk 

and the independence of the r.h.s from Xk (translational invariance) implies 

I (fx {Tkk{x))j = Lfc ^ In Z[J] . (B.5) 



(B.l) 

...On)j , 
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B.2 Translational invariance in 0-direction 



If we consider Eq. (1B.1|) for tiie interpolating field Tkk{vo), we choose (no summation 
over h) = e^Syk, and we integrate over the time-slices we obtain 

9^(^ok{xo)Tkk{yo))j = ai^{j(xo)([Tfefc(xo)+Z(xo)] Tfcfc(yo))j} 

+ tayo-xo)l(Tofc(yo))j . (B.6) 



By taking again Eq. (|B.ip for the field Tkkiuo) but for = eod^o, and by integrating 
over the full space we arrive at 



(B.7) 



do{Tkk{yo))j = - J dxoJ{xo)dQ{Tok{xo)Tkk{yo))j 
If we now put together Eqs. (|B.6p and (jB.7p we obtain 



d^o (Tkkiyo)). 



dxoJ{xo) Oq J(xo)((rfcfc(xo) + C{xo))Tkk{yo))j 



dpiyo) {Tok{yo))j . 



(Bi 



We are interested in applying the operator Pj^i_j,2n] defined in Eq. (j6.3p on both sides 
of the last equation. Since the Xq are all different from each other, the first term on the 
r.h.s. of Eq. ()B.8p does not contribute, and we obtain 



^[^..x-] {^o' {Tkk{yo))j + [dp{yo)] {Tok{yo))j} = 
where at most one of the Xq can coincides with i/q. In this case 

{j{yo)d^{Tok{yo))j} = 
and, since the derivative commute with Pr^,! 2.2^1, 



dK 



'0 {{Tkk{yo))j + J{yo) {Tok{yo))j} = 



(B.9) 



(B.IO) 



(B.ll) 



It is reassuring to notice that this equation can be obtained directly from Eq. ([2 
without the need of introducing Sj. 



By putting together Eqs. (jB.SP and (jB.lip we finally obtain Eq. ()6.9p . i.e. 

P[,i..,2„] {Lo (Tfcfc(yo)>j - ^ l^^t-^] - j J{wa) ^Qk{wQ))j] = , (B.12) 
where yo in this equation is different from all Xq. 
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C The Ward Identity ( 15. 3p in the free theory 

In this appendix we check exphcitly Eq. ()5.3p for the free theory. Other Wis, for instance 
Eqs. (14. Sp and (j5.4p . can be verified analogously. In a finite 4-dimensional volume the 
propagator 

{^{x)m) = D{x - y) (C.l) 



is given by 



D{x) 



V ^-^ rp' + 

V 



(C.2) 



where "p runs over a 4-dimensional momentum-space lattice = 2ti n^/L^ with /i 
1, ... ,4, and G ^. It satisfies the equation of motion 



[U-rr?] D{x-y) = -5^^'^ {x - y) 
where the periodic function 5^^"^ is defined as 



(C.3) 



(C.4) 



Thanks to the Wick theorem and to the invariance under parity 

(T^,(x)Tofe(2/)) = [d^d^Dix - y)] [d,dkD{x - y)] + [d^dkD{x - y)] [d,doD{x - y)] 

- V{ [dpd^Dix - y)] [d,dkD{x - y)] + n? [d^Dix - y)] [duD{x - y)] } , (C.5) 



where the derivatives are with respect to x. By using the equation of motion ()C.3p . it 
is easy to show that 

9^(r^,(a;)rofe(y)) = - [5o<5(P)(x-y)] \^d,dkD{x-y)\-[dk6^-^\x - y)] [d,doD{x-y)\ , 

(C.6) 

and by remembering that 



K^-Qk{x) = xoT^kix) - XkT^o{x), 
the symmetry of T^y implies (yo 7^ and yk 7^ 0, no summation over A;) 



(fxd^{K^,ok{x)TQk{y))=diD{z) 
By writing the r.h.s. as 



2=0 



dm 



z=Q 



(Too) - {Tkk) = dlD{z) 



2=0 



dlD(z) 



z=0 



then Eq. (jC.SP corresponds to Eq. (jS.Sp in the free theory, i.e. 



cfx d^{K^.ok{x)Tokiy)) = (Too) - {Tkk) 



2/0 7^ , yk^O ■ 



(C.7) 



(C.8) 



(C.9) 



(C.IO) 



18 



D Generating function for free and scale-invariant cases 



In this appendix we compute the generating function K{z) of the cumulants of the 
momentum distribution in the free theory, and check exphcitly the recursion relation 
()6.13p . We also compute the same quantity on the lattice to assess the magnitude of 
discretization effects in realistic computations. 

In a finite volume the shifted partition function is given by 



Z(z) 



HE 

p 



exp 



(-Lo^^p + «P • z) 



(D.l) 



where p runs over the 3-dimensional momentum-space lattice pi = Inmi/L with i=l,2,3 
and rrij G Z. As usual Up is the occupation number of the single-particle state labeled 
by p, see Ref. [1]. By performing the geometric sums, the cumulant generator is given 

by 

Ki-) = E'n . (D.2) 

p 



In the infinite volume limit 



E 



V 



(27r)3 



and rotational symmetry is restored. In the massless limit Up = [p| and, by integrating 
over the norm of p first, we obtain 



Kjz) 
V 



1 



1 



(1+^2^2)2 



(D.3) 



where the entropy density is given by s = 27r^r-^/45. In the massive case Wn = 
yjj^ + rr?, and we proceed by Taylor-expanding the logarithm and by performing the 
integration over the angular variables 



Kjz) 
V 



E 

n>l 



27r2n 



p^dp e 



-nLoujp 



sin(np|z|) 
np\z\ 



(D.4) 



By also Taylor-expanding the sinus, and by using the integral representation of the 
modified Bessel functions K^{x) we obtain 



V 



T(n + 



,5/2' E -2 E ^^TT)] (-2^'"^^^')" ^"+2(jWr) . (D.5) 



n> 



It is straightforward to verify that the cumulants generated from K{z) in Eqs. ()D.3p and 
(|D.5p satisfy the recursive relation (I6.13|] . For the massive case, the latter is implied by 



the recursive relation 



d Ky{ x) _ Ky+i{x) 
dx x^ 



among the modified Bessel functions. 
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D.l Free theory on the lattice 



Analogously to the continuum, the cumulant generator of the scalar free theory dis- 
cretized on the hypercubic in the standard way i.e. with 9^ replaced by the forward 
finite-difference operator, is 

^ ' ^ 1 - e-'f'O'^P 2 ^ cosh Lo - 1 

where 



2 asinh (\a^v^ + m^\ , p' = ^ E (^) ' ^^'^^ 



3 

„ . , / i /TTT- ^\ .9 4 , 

\'2 • ' I ' ' (7/ 

1=1 

By expanding in the lattice spacing a, in the thermodynamic limit it is easy to show 
that discretization effects are well-behaved and scale as expected, i.e. proportionally to 
(a/Lo)^. In practice one can use auj-p = 21nMp, Up = ^ay^p^ -|- m?+^^l + a?{-p^ + m?)/A, 

so that cosh(Lou;p) = ^(ttp^"^'^ -|- Up^^"^"), and carry out the sum over momenta nu- 
merically. 



D.2 Scale-invariant theories 

The functional form in Eq. (jP.SP is more generally valid than just for the free massless 
theory. The combination of scale invariance and of the recursive relation (j6.13p fully 
fixes the cumulant generator K{z). To show this we notice that scale invariance implies 

^{0,0,2n} = C{0,0,2n} T^^^^ . (D.8) 

It is then not difficult to solve the recursion relation (j6.13p for C{o,o,2n} to fiiid 

%o,2n} = (n + l)(2n)!|r2-. 
For z = {0, 0, z} the cumulant generator can thus be written as 

K{{0, 0, 4) = ^ E Vo,2n} , (D.IO) 

n=l ^ 

and by re-summing the series and using rotational invariance Eq. ()D.3P is obtained. A 
priori, one might have expected that in a scale-invariant theory, K{z) can be an arbitrary 
function of T\z\. However, the Ward identities of energy and momentum conservation 
fix its functional form unambiguously. 
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